Critical behavior in ultra-strong-coupled oscillators 
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We investigate the strong couphng regime of a linear x-x coupled harmonic oscillator system, by 
performing a direct diagonalization of the hamiltonian. It is shown that the x-x coupled hamiltonian 
can be equivalently described by a Mach-Zehnder-type interferometer with a quadratic unitary 
operation in each of its arms. We show a sharp transition of the unitary operation from an elliptical 
phase rotator to an elliptical squeezer as the coupling gets stronger, which leads to the continuous 
generation of entanglement, even for a significantly thermal state, in the ultra-strong coupled regime. 
It is also shown that this critical regime cannot be achieved by a classical Hookian coupling. Finally, 
the effect of a finite-temperature environment is analyzed, showing that entanglement can still be 
generated from a thermal state in the ultra-strong coupled regime, but is destroyed rapidly. 

PACS numbers: 42.50Ct, 42.50.St, 03.67.Bg 



I. INTRODUCTION 

Strong coupling, which is one of the key ingredients for 
the manipulation and control of quantum systems, makes 
quantum-mechanical predictions distinct from classical 
ones. In fact, the definition of strong coupling in quan- 
tum mechanics has evolved over the years. In atomic 
and optical physics, strong coupling was defined to show 
Autler-Townes splitting [T] where atom-field interaction 
is stronger than the atomic decay rate. The next age of 
strong coupling was defined in cavity quantum electro- 
dynamics (QED), where an atom- field interaction strong 
in comparison to the cavity and atomic decay rates can 
be realized. A series of important experiments have been 
performed to show the nonclassical nature of the atom- 
field interaction for Rydberg atoms in microwave cavities 
[H [3] . This regime was recently also attained with a cav- 
ity optomechanical system [J]. But with a system con- 
sisting of superconducting Josephson junctions coupled 
to microwave fields in stripline resonators, the interac- 
tion strength can be made even stronger than the atom- 
field coupling realized previously in cavity QED [51 E]. 
These so-called circuit QED systems can reach coupling 
strengths comparable to the transition frequency of the 
relevant qubit system. In these systems, we find new 
phenomena such as the breakdown of the rotating-wave 
approximation and the emergence of a unique deep strong 
coupling regime 7]. 

A system of coupled oscillators has been of interest in 
various contexts, as several physical systems are repre- 
sented by harmonic oscillators such as nanomechanical 
oscillators, electromagnetic fields and a bunch of two- 
level systems. Hopfield studied the quantization of an 
electromagnetic field in a dispersive medium using a cou- 
pled harmonic oscillator model . Recently, such a sys- 



tem of coupled oscillators has been considered as a pos- 
sible mechanism to generate entanglement in an array 
of nanoelectromechanical devices [Qllll). In this paper, 
we study ultra-strong coupling between harmonic oscilla- 
tors where the coupling rate g is comparable to or larger 
than the natural frequencies loj of the harmonic oscil- 
lators. Strongly coupled oscillators, where the counter- 
rotating terms are not negligible have been studied in 
the context of twin-photon generation |12j . The study 
of such systems are supported by recent experimental 
progress wherein such a regime is becoming accessible, 
for instance in the interaction of multiple quantum well 
structures with terahertz electromagnetic fields [T^. In 
this paper, we go even further than the strong coupling 
regime, and find a critical point in the coupling strength 
at which the evolution of the oscillator states radically 
changes. An equivalent regime was studied for the inter- 
action between an atom and a field [7l [M] , where pho- 
ton number wavepackets were found to experience col- 
lapse and revivals across parity chains in Hilbert space. 
The model that we consider here, is quite different in the 
sense that we have a bosonic system living in an infi- 
nite dimensional Hilbert space, so that there is no such 
parity-defined dynamics. However, the model is exactly 
solvable and we find critical effects from the dynamics 
of the coupled oscillators, which get manifested in their 
entanglement behaviour. 

We consider a system of two bosonic oscillators inter- 
acting via a quadratic interaction, whose hamiltonian is 
described by 
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where the dimensionless quadrature operators satisfy the 



canonical commutation relation, i.e. 
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instance, the Dicke model of superradiance [IS] is ap- 
proximately described by such a hamiltonian. In the con- 
text of cavity optomechanical systems, such a linearized 
hamiltonian describes the dynamics when the mechanics 



is being driven by a strong coherent field [T5]. In fact, 
this Unearized model has been investigated previously in 
the optomechanics context [17] in the rotating wave ap- 
proximation (RWA), which is possible when the system 
is relatively weakly coupled. 

In this paper, we will show that in the ultra-strong 
coupling regime, the non-RWA terms turn out to be not 
only important but also to critically change the oscillator 
interaction. After identifying a unitary transformation 
exactly diagonalizing the hamiltonian in (fTj) , we show its 
equivalence to a Mach-Zehnder-type setup with unitary 
operations in its arms, which explains the generation of 
entanglement in the model. At the critical point when 
the coupling becomes ultra-strong, we observe the tran- 
sition of the unitary operation from an elliptic rotator to 
a squeezer. This transition is clearly reflected in the en- 
tanglement dynamics of the oscillators. It is interesting 
to note that the classical Hookian model for classical har- 
monic oscillators is not able to reach this critical point. 



II. DIAGONALIZATION OF THE 
HAMILTONIAN 



The diagonal hamiltonian is then given by. 
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where we denote c — cos y/AB and s = sin y/AB. Each 
normal mode is described by a quadratic hamiltonian, 
which looks very similar to a standard free hamiltonian 
for an oscillator. But using the condition Q, it is found 
that the last term in (pi) is zero (negative) when g = gc 
{g > gc), so that these are standard harmonic oscillators 
only when g < gc- The first (second) mode in ^ is 
associated with E+ {E-), which we call the '-I-' (^^— ') 
mode. When g = ga the '— ' mode does not have a 
bound spectrum, but is rather a free particle, while the 
'+' mode is still a harmonic oscillator. Increasing the 
coupling strength further, for g > g^ the anomalous '— ' 
mode is dynamically unstable since it is driven by a force 
derived from the inverted harmonic potential. So, the 
earlier observation of £'_ going imaginary is reflected by 
this qualitative change in the dynamical behaviour of the 
system, as the coupling strength crosses its critical value 



One can codify the operator hamiltonian (IT]) in terms 
of a hamiltonian matrix. 
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expressed in the basis defined by the column vector of 
quadratures R = {xi,X2,Pi,P2y , so that H = ^WHR 
where we also define the diagonal block matrices corre- 
sponding to positions and momenta. 

The usual diagonalization of such a quadratic hamilto- 
nian proceeds by defining normal mode operators J7g [q = 
±) and their hermitian conjugates, which are linear com- 
binations of the bare quadratures satisfying [Clq , H\ = 
Eqflq. Specifically, the normal mode energies are found 
to be E'^ = wj -I- wf ± \/{ujI + uj\Y + 4a;icj2(5^ — W1W2), 
from which we already observe an interesting regime 
where E- is imaginary, for g > gc^ ^JuJ\uj2. 

The transformation to normal modes is just a sym- 
plectomorphism [18 . The hamiltonian matrix can be 
diagonalized by a symplectic matrix S £ Sp{A, M) such 
that it is block diagonal of the form S — Sx (B Sp. The 
diagonal hamiltonian matrix is then H' — S^HS — 
(SJHxSx) ® (SjHpSp). From the general conditions to 
be satisfied by symplectic transformations, we get that 
for Sx,Spe SL{2,R), Sj = S'^ The unitary represen- 
tation of S in the Hilbert space of the system is, 



T = exp [i{AxiP2 - BX2P1) 
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parameterized by real constants A and B. This diago- 
nalizes H to H' = THT^ iff. 
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III. INEQUIVALENCE TO CLASSICAL 
COUPLED OSCILLATORS 

An obvious intuition one would have from thinking 
about the classical oscillating systems is that the model 
(fTl) must be the quantized version of the classical Hookian 
hamiltonian for coupled oscillators. 
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where m is the common mass of the particle, w the com- 
mon natural frequency and Q the Hookian coupling rate. 
The oscillators are assumed to be identical for simplic- 
ity (a generalization is straightforward). Introducing di- 
mensionless quadratures, a;i_2 = ("^^^0)^^^1,21 Pi, 2 = 
(ttiojo)"^ ^^1,2 with the renormalized frequency wq — 
\/ uj^ + t/2^ the hamiltonian becomes. 
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where G 
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is the renormalized coupling strength. 



So, it is clear that in the quantum version of the classical 
Hookian problem, the ultra strong coupling regime (here, 
G > 1) is not possible to achieve. 

The ultra strong coupling transition of the form ex- 
hibited by the hamiltonian ( ll) seems to be a property of 
such quantum coupling models, and not something that 
is seen easily in the standard harmonic oscillator model. 
In particular though, for optomechanical problems, the 
linearization of the actual radiation pressure interaction 



{p1 + x1)x2 (which leads to i;ii;2) is valid only for a short 
interaction time. Otherwise, the displacement of the me- 
chanics gets too large and the approximation leading to 
the interaction hamiltonian breaks down in the strong 
coupling regime g > Qc 



IV. ENTANGLEMENT DYNAMICS 

In the resonant case, i.e. Wi = uj2 and bo A ~ 
B = J, the transformation operator T in Eq.(^3| be- 
comes exactly the same as the 50:50 beam splitter op- 
erator. Thus the evolution operator U{t) = e^'-^* 
corresponding to the hamiltonian (fTl) can be factor- 

ized into U{t) = f^e-'"-*e-''"'+^f, where H'^ = 
^ {'-^1.2'Pi 2 + {9c i g)x1 2) ■ It is now obvious that in the 
resonant case, the dynamics is equivalently represented 
by a Mach-Zehnder interferometer, with 50:50 beam- 
splitters as shown in Fig. [T] but with some active op- 
erations in both arms [25]. 

To determine the kind of active operations, consider 
the individual evolutions generated by H± which are of 
the quadratic form, Hg = ^(a^a;^ -I- /3qPg) {q = ±). The 
evolution of the quadrature operators induced by such a 
hamiltonian is 
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The transformation shows the behavior of an elliptical ro- 
tator where the product of the quadrature scaling factors 
is unity. It preserves the quadrature uncertainties, lead- 
ing to at most, squeezing of the mode. Importantly, if ei- 
ther Uq or /3g is negative, the elliptical rotator becomes an 
elliptical squeezer ~ this is precisely the situation for the 
'— ' mode when the coupling is ultra strong, i.e. g > ge- 
lt is also worth noting that in the general non-resonant 
case, both the beam-splitting and two-mode squeezing 
terms appear in T, meaning that instead of the passive 
beam-splitters at the entry and exit ports, one will have 
another active device in the equivalent picture. 
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FIG. 1. The evolution due to hamiltonian ([I| is equivalent 
a Mach-Zehnder setup for the resonant case, such that one 
arm of the interferometer experiences elliptic rotation and 
the other arm undergoes elUptic rotation (squeezing) ii g < Qc 
(g > gc). BS; Beam spUtter, M: Mirror. 



From the equivalent Mach-Zehnder picture, the source 
of entanglement is obvious - squeezing in the input fields 
is a necessary condition for Gaussian entangled output 
from a beam splitter [57. The interaction leads to sig- 
nificantly more entanglement when g — gc, than the case 
where g < gc- To explicitly study dynamic entangle- 
ment, we consider each harmonic oscillator initially in its 
thermal equilibrium so that the total density operator 
p(0) = pJ^-Sip^^, where pj'^ = ^ exp(-^(x|-fp2)) with 
the canonical partition function Zj and mean excitation 
number rjj = (e ^ — l)"'^. Sj = uij/ksTj at temperature 
Tj with the Boltzmann constant ks- 




FIG. 2. (Color online) (a) Log negativity, and (b) the seralian, 
for the thermal state with rj^ — 0,rj2 — 1 for various coupling 
strengths, at cJi = 5uJ2- 




FIG. 3. (Color online) Dynamic entanglement at g — g^ for 
increasing degree of thermality when one oscillator is initially 
(a) in a pure state {rji = 0) and (b) in a thermally excited 
state {rji = 5). The two oscillators are resonant. 

Fig. p[a) shows the dynamic entanglement, measured 
by the logarithmic negativity [53J [53] , for a thermal state 
with a single average excitation (see the appendix |A] for 
explicit formulas of the log-negativity) . We observe that 
for g < gc, the entanglement has an oscillatory behaviour, 
with typical periods related to the frequencies of the nor- 
mal modes. When g ~ g^ the entanglement monotoni- 
cally increases, disregarding the small oscillations due to 
the '-I-' normal mode. The same behaviour is observed 
for larger values of the coupling constant (g > gc) in 
Fig. I2jb), where we plot the "seralian" A [26] which, in 
the unitary case, is a sufficient entanglement monotone 
(see Appendix IA] for details) . The entanglement keeps 
increasing with time (neglecting the small oscillations). 



and in particular at a fixed time, the larger is the value 
of the coupling constant g, the larger is the entanglement 
achieved. This particular behaviour can be explained also 
by looking at the eigenvalues of the Hamiltonian matrix 
H. For g < gc, the matrix is positive definite and, as 
proved in p7j . the corresponding unitary operator will 
recur to the identity for a given time; for larger values 
of g, the matrix H is not positive definite anymore, the 
recurrence property of the evolution is lost and the en- 
tanglement increases with time. 

Even when the initial state is significantly thermal, if 
the system is critically coupled, entanglement can be gen- 
erated, as depicted in Fig. [3j In fact, we also observe 
that the purity of one oscillator is an important factor 
in the generation of entanglement. When one oscillator 
is pure, the maximum degree of entanglement achieved 
remains about the same regardless of the temperature of 
the other oscillator. The importance of single-system pu- 
rity was also observed in some other interaction models 




g=0.8g,;n=(l,10) 

g=g,;n=(l,10) 

g=1.5g,;n=(l,10) 

-— g=0.8g,;n=(l,5) 
.n=(l,5) 



-— g=1.5g,;n=(l,5) 



FIG. 4. (Color online) Dynamic entanglement for the dis- 
sipative case, for each of the three regimes. Here the oscil- 
lator 1 (2) is coupled to an environment of ni = 1 at a rate 
71 = 0.01a;i (72 = 0.25a;2). 



A. Dissipative dynamics 

Let us consider now the case where the system in- 
teracts with a noisy environment. Because of the en- 
vironment interaction we observe that also super-critical 
regime g > gc can stop to be dynamically unstable. As- 
suming that each bare oscillator interacts with its respec- 
tive thermal environment with mean excitation n^, un- 
der the Born-Markov approximation, the evolution of the 
system density matrix is governed by the Kossakowski- 
Lindblad master equation. 



dp 
~di 






(7) 

where C[d]p ^ 20 pO'^ - O'^ 6 p ~ pO'^ 6 is the Liouvillian, 
while 7j is the rate at which the j*'^ mode is coupled 
to the environment (see Appendix p] for the analytical 
solution of I?])) 

The logarithmic negativity Ejq can be evaluated, and 
the results are plotted in Fig. |4]for various values of the 
coupling strenght. We observe that, as expected, also 
in the super-critical regime g > gc, the entanglement 
reaches a maximum value and then starts to decrease 
and eventually reaches zero for a given time to. Not sur- 
prisingly, larger values of the coupling constant and lower 
values of the noisy parameters rjj , jj and fij , correspond 
to larger maximum values of entanglement and to larger 
values of the zero-entanglement time to. 



V. CONCLUSION 

We have investigated the ultra-strong-coupled regime 
of a linear x — x coupled quantum harmonic oscillator 



system. Three different regimes of qualitatively different 
dynamical behaviour can be identified from the normal 
modes of the system. In particular, in the super-critical 
regime g > gc, the '— ' mode has an anomalous acceler- 
ating motion, which can be understood as a consequence 
of a lack of translation symmetry in the hamiltonian (II]) , 
so that momentum conservation is not explicitly guar- 
anteed. In the classical Hookian hamiltonian ^ such a 
symmetry is manifest, so that no dynamical anomaly is 
allowed in the strong coupling regime. 

We further find that these three regimes characterized 
by their dynamical behaviour, is also uniquely charac- 
terized by the dynamics of entanglement; in particular, 
that the super-critical regime is one where entanglement 
can be unboundedly generated from even a highly ther- 
mal state. A similar behavior was studied in the con- 
text of the superradiant phase transition of the Dicke 
model, whose effective hamiltonian is of the form (fTl) af- 
ter a Holstein-Primakoff transformation of a set of SU(2) 
operators into a bosonic operator. 

Finally, we also note that the super-critical regime 
exhibits a finite entanglement even though the system 
maybe coupled to a decohering Markovian environment 
at a non-zero temperature. 
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Appendix A: Log-negativity for Gaussian states 

Given a two-mode Gaussian state, its entanglement 
properties are fully characterized by its covariance matrix 
(T, which can be written in terms of 2 x 2 matrices as 



7^ (^2 



(Al) 



The log- negativity [21| of the two-mode state can be eval- 
uated as Em{q) = max[0, — log(2i'_)], where U- is the 
lowest symplectic eigenvalue of the corresponding partial 
transposed state. In formula we have 



vl = 



A-^A2-4Det[(T] 



where the "seralian" A is defined as 



A = Det[cri] + Dot [0-2] - 2Det[7] 



(A2) 
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and where Det[A] denotes the determinant of the matrix 
A. It is known that the logarithmic negativity Ejsi is an 
entanglement monotone |1S]. One should also note that 
Det[(T] is directly related to the purity of the state, since 
Tr[£)2] — l/{2^JT)et[(T]). As a consequence, unitary evo- 
lution leaves Det[<T] invariant, so that the log-negativity 
is monotonous with A. For this reason, when we consider 
the entanglement dynamics due to non-dissipative evolu- 
tion, we may take the seralian A to be an entanglement 
monotone. 



Appendix B: Solving the master equation 

The master equation (l7| can be cast into a c-number 
partial differential equation (PDE) by introducing the 
symmetric ordered characteristic function 

X{ai,a2]t) = Tr[I)i(ai)I)2(a2)/5(i)], 

where Dj{aj) = exp{ajaJ- —a*aj). We choose to express 
the PDE in terms of the quadratures parameters Xj,pj 
defined via aj = 4= {xj + ipj ) , so that we get the Fokker- 
Planck equation. 
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where we define the vector r 
sociated gradient operator dr 
following matrices: 
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and where H here is a permutation on the hamiltonian 
matrix (pi) corresponding to the difference in the order 
chosen for the quadrature variables. The arrows on the 
gradient operators denote the direction in which the dif- 
ferential operators act. 

As shown above, the system Hamiltonian only affects 
elliptic rotations or elliptic squeezing in the course of the 
dynamics, so that an initial Gaussian state localized at 
the origin in phase space remains so always; and since the 
bath is in thermal equilibrium, it too does not affect any 
finite displacements in phase space, but only scale and 
rotation changes. Thus, for our relevant case, without 
loss of generality, we choose the ansatz. 



X(r;t) =exp ( --r"^fT(f)r 



where now the time dependence is carried in the covari- 
ance matrix. Substituting this into (Bl ), and then iden- 



tifying the coefficients of the various bilinear products, 
we get the equation of motion for the covariance matrix 
viz., 
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This has the solution. 



a{t) = K{t) { cr(O) + / K{-t)TK ' (-r) dT\K\t), 



(Bl) where K{t) = exp {tH - \t\ 
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